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ABSTRACT 



Aims. Procyon A, a bright F5 IV-V Sun-like star, is justifiably regarded as a prime asteroseismological target. This star was repeatedly 
observed by MOST, a specialized microsatellite providing long-term, non-interrupted broadband photometry of bright targets. So far, 
the widely anticipated p modes eluded direct photometric detection, though numerous independent approaches hinted for the presence 
of signals in the f~ 0.5 - 1.5 mHz range. 

Methods. Implementation of an alternative approach in data processing, as well as combination of the MOST data from 2004 and 
2005 (264 189 measurements in total) helps to reduce the instrumental noise affecting previous reductions, bringing the 3tr detection 
limit down to ~5.5 part-per-million in the / = 0.8 - 1.2 mHz range. 

Results. This enables to cross-identifiy 16 p-mode frequencies (though not their degrees) which were previously detected via high- 
precision radial velocity measurements, and provides an estimate of the large spacing, Sv = 0.0540 mHz at / ~ 1 mHz. The relatively 
low average amplitude of the detected modes, a = 5.8 ± 0.6 ppm, closely matches the amplitudes inferred from the ground-based 
spectroscopy and upper limits projected from WIRE photometry. This also explains why such low-amplitude signals eluded the 
direct-detection approach which exclusively relied on the MOST 2004 (or 2005) data processed by a standard pipeline. 
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1. Introduction 

Procyon A (F5 IV-V) has served as a prime target in many 
asteroseismological campaigns (see the reviews by Bedding 
& Kjeldsen, 2006,2007, and references therein). While 
high-precision measurements of radial velocities revealed 
the presence of variability caused by solar-like oscillations 
(Mar tic etal. 19 99) and provided a preliminary identifi- 
cation of p-mode frequencies (Eggenbe rgeret al. 2004| 
Marti c et al. 20041 ILeccia et al. 2007t , the photometric studies 
met either limited success dBruntt et al. 20051 ) or outright 
null-detection dMatthews et al. 20041 IGuenther et al. 2 007). 
The null-detection result was hotly debated ever since 
(|Bedding et al. 2005| |Regulo & Roca Cortes 2005[). Attempting 



to resolve the controversy, we applied an alternative data- 
processing approach to the abundant photometric data acquired 
by MOST in 2004-2005. 



2. Observations and data processing 



The MOST satellite ( Walke retal. 20031 ) is a 15-cm telescope 
acquiring high-precision, non-stop CCD photometric data on 
bright objects through a broad-band, 350-700 nm, optical fil- 
ter. MOST observed Procyon A in 2004 (32-day non-stop ob- 
servations, 231 524 0.9 s exposures: see the results discussed in 
Matt hews et al. 2004] ), 2005 (taken over 17 contiguous days with 
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* MOST is a Canadian Space Agency mission, operated jointly by 
Dynacon, Inc., and the Universities of Toronto and British Columbia, 
with assistance from the University of Vienna. 



the general setup of the experiment following the 2004 cam- 
paign: Guenthe r et al. 2007] ) and 2007 (data are being processed: 
J. Matthews, 2007, priv. comm.). We took the publicly avail- 
able raw data from 2004 and 2005 and processed them using an 
alternative algorithm (see the description in lAerts et al. 2 006 ). 
While generally conforming to the approach of the standard 
processing pipeline applied to the 2004-2005 data, the alterna- 
tive algorithm treats the highly variable stray light, the main 
source of instrumental noise, in a more flexible manner. In 
addition to the iterative subtraction of scattered light, the al- 
gorithm accounted for less substantial sources of instrumen- 
tal noise, such as sensitivity drift caused by slow variations of 
the CCD temperature. The alternative approach resulted in a 
lower point-to-point scatter, s =440 part-per-million (ppm), vs. 
s =500 ppm ( Matthew s et al. 2004] ) (cf. the upper section of 
Fig 1). This should be compared to the projected level of the 
Poisson-dominated instrumental noise s,-„ =260 ppm. The alter- 
native approach also reduces the average (noise-dominated) am- 
plitudes and scatter in the frequency spectrum: e.g., for f=1.05- 
1.25 mHz one obtains a(f) = 2.8 ± 1.7 ppm in the older re- 
duction vs. a(f) = 2.6 +1.4 for the new algorithm (Fig. [l}. 
The modified approach, when efficiently combining the two 
available data streams (SDS1 and SDS2: see the definition in 
Aerts et al. 20061 ). proves to be even more effective with the 
2005 data where the better planning of experiment allowed more 
straightforward suppression of the stray-light component, pro- 
viding s ^315 ppm, which only moderately exceeds the ex- 
pected level of Si„ =250 ppm. Predictably, the lower noise 
level brings down average amplitudes in the frequency domain: 
«(/) = 1.4+1.6 ppm in the f= 1.05- 1.25 mHz range, twice lower 
than in the data from 2004. On the other hand, the standard 
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Fig. 1. The upper section shows the new spectrum (black lines) derived from the stand-alone 2004 data and the previously published 
(red lines) spectrum derived from the same data. The orbital harmonics and their n* 1/d sidelobes are marked by vertical solid-line 
and dashed-line arrows, respectively. Middle section: the Fourier amplitude spectrum calculated from the combined 2004-2005 data 
set processed with the alternative algorithm (black lines), compared to the published (Matt hews et al. 20041 red lines) spectrum as 
a reference. Blue arrows mark the identified p-modes. Bottom section: an artificial spectrum calculated from the shuffled MOST 
2004+2005 data with added stochastically-excited sinusoidal signals. Their frequencies are taken from Table 1. For simplicity, all 
the signals are assigned equal amplitudes a=10 ppm and lifetimes L=7 days (red lines), and a=15 ppm, L=3 days (black lines). Blue 
arrows mark the identified p-modes. In all sections the horizontal blue lines denote the 3cr detection limits. 



approach (Guentheret a l. 20071) . relying uniquely on the SDS2 
channel, provides a noise level comparable to the 2004 set. 

During the processing of the 2004 data we kept 
the time-binning approach similar to the original scheme 
(Matthews et al. 2004b : this provided a valuable point-to-point 
comparison between the two different approaches and retained 



214279 flux measurements. For the 2005 data we slightly 
changed the approach, by time-binning the SDS1 observations 
centered on the less frequent SDS2 measurements. On average, 
taking into account the higher rejection rates, the SDS1 data 
were ~ 5.5 times more abundant in the 2005 set. This time- 
binning approach resulted in 49 910 flux measurements. 
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3. Results 

Our first step is to compare the frequency spectra calculated by 
the two data-processing algorithms. As anticipated, the alterna- 
tive approach, while providing a frequency spectrum of an over- 
all similar appearance, reduces the level of instrumental noise 
across the spectrum, save a few rare exceptions related to the 
harmonics of the orbital period (black arrows in Fig. |T}. 

The next step aims at an optimal combination of the data. 
A host of potential problems stems from: (a) the different du- 
rations of the 2004 and 2005 runs; (b) different planning of 
the experiments resulting in different levels of stray light; (c) 
presumably unstable nature of the signals. In the time do- 
main we experimented with different lengths of segments de- 
rived from the original time series, varying them from 3 days 
to 2 weeks and producing a weighted combination of the fre- 
quency spectra derived from the segments. This allowes us to 
see that the typical lifetime of signals does not exceed ~ 1 
week. However, the segment-based approach comes at high 
price, raising in the region of interest, f=0. 5-2.0 mHz, the de- 
tectability limits to an unacceptably high level, a > 10 ppm. 
Hence, we finally relied on the simplest and most straightfor- 
ward approach, by combining the unweighted 2004 and 2005 
data into a single set (Fig. Q]). Obviously, this results in a fur- 
ther improvement of sensitivity: considering the 3-sigma level 
counted from the average amplitude of the noise-dominated sig- 
nal, one is able to detect coherent oscillations with a > 6.0 
ppm (f=0.7-0.8 mHz), a > 5.5 ppm (f=0.8-1.2 mHz), a > 5.0 
ppm (f =1.2-1. 6 mHz), a > 4.5 ppm (f>1.6 mHz), to be com- 
pared to a > 7.0 ppm (f =0.5-2.0 mHz: IMatthews et al. 20041 1 
and a > 10—11 ppm (Guen ther et al. 20071) . Beyond the over- 
all improvement of sensitivity, such straightforward combina- 
tion of the 2004 and 2005 data should substantially increase the 
chance of a positive detection which depends on the presum- 
ably short lifetimes of the coherent signals (p-modes) in Procyon 
A dLeccia et al. 20071) and the stochastic nature of the mode- 



excitation events (Kjeldsen & Bedding 1995 1 



As a first test of presence of any periodic variations, we 
derived a Fourier amplitude spectrum (AS) of the combined 
2004+2005 data in the f=0.5-2.0 mHz domain specifically tar- 
geted in numerous observing campaigns. Then we calculated 
the number of peaks exceeding the 3cr detectability levels (see 
above). We also produced 10 samples via random shuffling of 
the real data and calculated the incidence of significant signals in 
the artificial ASa. We normalized the real data by the outcomes 
of simulations and show them in Fig. [2] There is a surplus of 
3<r detections around f~0. 8-1.4 mHz, which closely matches the 
spread and location of the power excess determined from nu- 
merous independent spectroscopic runs (see ILeccia e t al. 2007 
and references therein). 

Further testing brings even more encouraging results. We 
select a prominent feature at f= 1.072 mHz (Table Q] and 
Fig. [TJ and search for a comb-like pattern (Kjel dsen et al. 1995) 
in the immediate vicinity of the peak. Scanning the ad- 
jacent region with Sf - 0.0001 mHz steps, we find 
a maximum response at f= 1.0726 mHz (Fig. |3), related 
to the large frequency separation of Af = 0.0540 + 
0.0001 mHz, in good agreement with Af = 0.0536 - 
0.0559 mHz consistently provided by all previous spectro- 
scopic campaigns (Eggenberger et al. 2004 IMartic et al. 2 004; 
ILeccia et al. 2007t . The clear presence of a regularly-spaced 
signal in the MOST-2004 data was noted on multiple occa- 
sions: A/ = 0.0545 mHz ( |Regulo & Roca Cortes 2005) , Af ~ 
0.0548 mHz (J.Matthews, priv. comm. 2005), also by the ref- 
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Fig. 2. The numbers (iV,) of significant (at 3<x level) peaks in 
the 2004+2005 AS normalized by the corresponding numbers 
(N s ) derived from artificial data sets. The vertical bars show the 
±lcr errors of the ratios. The horizontal dashed line provides a 
reference level of 1 (the case when all peaks are generated by 
noise). 



eree of the article. Can the regular spacing of signals be caused 
by an unaccounted component of the stray light contributing to 
the n * v or t,/3 orbital harmonics with 0.0548 mHz separation? 
Apparently not, as the newly derived Af is more than sufficiently 
distanced from the instrumental component. Moreover, we find 
only a very weak trace of the instrumental signal in the comb- 
generated response at f= 1.0726 mHz and its complete lack at 
f=1.6424 mHz (Fig. [3] note that f= 1.6424 mHz matches an or- 
bital harmonic). 
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Fig. 3. The comb-response functions obtained for the central 
frequencies f= 1.0726 mHz (full line), f=0.9158 mHz (dashed 
line) and f= 1.6424 mHz (dashed line, offset for clarity). The 
tallest peak corresponds to a large separation Af = 0.0540 ± 
0.0001 mHz, with the error derived from the FWHM of the peak. 



The presence of a regular 'comb' of periodic signals around 
f= 1.072 mHz prompted the final step of the analysis. We 
compiled a list of frequencies presumably related to p-modes 
(Egg enberger et al. 2004| IMartic et al."2004l ILeccia et al."2"0"07j l 
and searched for corresponding significant peaks in our 
2004+2005 AS. We limited our search to the immediate sur- 
roundings of the published frequencies, i.e., to the +5/^0.002 
mHz intervals, where Sf is an uncertainty related to frequency 
resolution in the AS derived on ~1 week-long runs. Note that 
a formal frequency resolution in the combined 2004+2005 AS 
reaches 3 * 10~ 5 mHz. We recorded the maxima which were 
very close to or above the corresponding 3<x levels for the 
given frequencies. Facing the rapidly rising level of instrumen- 
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Table 1. Detected p-modes 



f [mHz] 




Amplitude, 


a 


[ppm] 


Refei 




2004+2005 


2004 


2005 


2004 1 




0.7041 


7.1 


7.4 


4.7 


8.0 


3 


0.7206 


6.0 


6.0 


5.6 


_ 


2,4 


0.7525 


(5.6) 


6.2 


5.3 


(6.6) 


3 


(0.7692) 


6.2 


6.3 


4.0 


7.5 


3 


0.8027 


(5.0) 


5.7 


4.0 


_ 


2.3 


0.9158 


6.1 


6.6 


5.3 


9.0 


2,4 


1.0232 


5.8 


6.0 


3.6 


_ 


2,3,4 


(1.0726) 


6.2 


9.2 


8.1 


7.6 


3 


1.1136 


5.4 


5.7 


3.9 




4 


1.1837 


5.6 


6.5 


2.1 




2,3 


1.1959 


6.9 


7.8 


3.5 


(6.2) 


2,4 


1.2158 


5.0 


5.4 


4.2 




3 


1.2415 


6.5 


7.6 


5.2 


(5.2) 


3 


1.2692 


5.9 


6.1 


3.3 


(5.7) 


2,3,4 


1.3547 


5.2 


5.9 


3.3 


(6.0) 


3 


1.5590 


4.9 


4.9 


4.4 




2,4 



them in Fig. 1 . In both shown cases, a an — 1 ppm, L-l days, and 
a art -\5 ppm, L-3 days, one may detect up to ~ 80 - 90% of the 
artificial signals. The a an = 10±ppm case provides the average 
amplitude < a an >- 5.0 ±0.8 ppm, and a ar ,=15 ppm case results 
in < a ar , >- 7.1 ± 1.2 ppm and a slightly better detection rate, 
thus proving the starting assumption about the short lifetime of 
the modes. On the other hand, the relative low signal/noise ratio 
of the data and, presumably, interaction of the stochastically ex- 
cited modes affects the derived frequencies of artificial signals 
(fart)- < [fobs - fartl >~ 9 * 10~ 5 mHz, to be compared to the 
formal frequency resolution 6f — 3* 10~ 5 mHz in the combined 
2004+2005 data.. 
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operated by the Hertzberg Institute of Astrophysics, national Research Council 
of Canada). Numerous suggestions of the anonymous referee helped to improve 
presentation of the main results. 



References 

Aerts, C, et al. 2006, ApJ, 642, 470 
Bedding, T.R., et al. 2005, A&A, 432, L43 

Bedding, T.R., & Kjeldsen, H. 2006, ESA Special Publication, 624, 25 
Bedding, T.R., & Kjeldsen, H. 2007, Comm. in Asteroseismology, in press 

(arXiv:0705.1376l 
Bruntt, H., et al. 2005, ApJ, 633, 440 

Guenther, D.B., et al. 2007, Comm. in Asteroseismology, 151, 5 

Eggenberger, P., et al. 2004, A&A, 422, 247 

Kjeldsen, H., & Bedding, T.R. 1995, A&A, 293, 87 

Kjeldsen, H., et al. 1995, AJ, 109, 1313 

Leccia, S., et al. 2007, A&A, 464, 1059 

Martic, M, et al. 1999, A&A, 351, 993 

Martic, M., et al. 2004, A&A, 418, 295 

Matthews, J.M., et al. 2004, Nature, 430, 51, Erratum: 430, 921 
Regulo, C., & Roca Cortes, T. 2005, A&A, 444, L5 
Walker, G., et al. 2003, PASP, 115, 1023 



1 IMatthews et al. 20041 

2 Eggenberg er et al. 2004| 

3 Ma rtic et al. 20041 

4 ILeccia et al. 20071 raw data and fits. 



tal+intrinsic (granulation?) noise, we disregarded all positive de- 
tections with frequencies below f=0.7 mHz. We also eliminated 
all positive detections (4 cases) which match high harmonics of 
the orbital period and their 1/d side lobes. We present the re- 
sults in Table [TJ where in parentheses we retain some promi- 
nent signals despite their relatively large deviations from the 
published frequencies. For a reference, we also provide the am- 
plitudes of signals in the separate subsets from 2004 and 2005, 
even if they do not exceed detection limits. We compare our pos- 
itive 2004+2005 detections with signals exceeding (or being rea- 
sonably close to) the detectability limits in the previously pub- 
lished MOST-2004 data (IMatthews et al. 2004K finding ~ 55% 
of positive matches, both in frequency and amplitude. On aver- 
age, one may expect 2.4 + 1.2 (~ 20 in total) noise-generated 
signals exceeding the imposed 3<x threshold in an 0. 1 mHz-wide 
bin placed in the region of interest, f=0.7-1.6 mHz. Hence, the 
chance for an accidental match between a noise-generated peak 
and a pre-determined signal registered in the previous spectro- 
scopic campaigns ( Eggenberger et al. 2004 Martic et al. 2004) 
is < 45%. I.e., at least 50% of the frequencies provided in TableQ] 
should be geniune. The 16 signals identified as p-modes, plus 20 
noise-generated peaks do not account for all positive detections, 
46 signals in the f=0.7-1.6 mHz range. The origin of the 10 re- 
maining peaks is unclear. 

The average amplitude of the recovered signals, a - 5.8 ± 
0.6 ppm, closely matches the a = 7.3 + 1.1 ppm level in- 
ferred from the ground-based spectroscopy ( ILeccia et al. 2007l l 
and a = 8.5 ± 2 ppm upper limit from WIRE photometry 
(Bru ntt et al. 20 05). Such low-amplitude signals were inevitably 
swamped by instrumental noise in the previous reduction of the 
stand-alone data from 2004 and 2005. The unusually low level of 
the detected signals may be related to the anticipated short life- 
time of the modes, 2.0 + 0.4 d ( ILeccia et al. 20071 . Indeed, one 
may re-shuffle the original 2004-2005 MOST data and add 16 
sinusoidal signals with the frequencies matching the identified 
p-modes from Table 1 . For simplicity, we assign to the signals 
the same amplitude a art and lifetime L, let the phases of the sig- 
nals vary randomly, then calculate the resulting ASa and show 



